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ABSTRACT 

The Planetary Nebula (PN) population of M33 is studied via multi-fiber spectroscopy with 
Hectospec at the MMT. In this paper we present the spectra of 102 PNe, whereas plasma diag- 
nostic and chemical abundances were performed on the 93 PNe where the necessary diagnostic 
lines were measured. 

About 20% of the PNe are compatible with being Type I; the rest of the sample is the progeny 
of an old disk stellar population, with main sequence masses M<3M0 and ages t>0.3 Gyr. 

By studying the elemental abundances of the PNe in the M33 disk we were able to infer 
that: (1) there is a tight correlation between 0/H and Ne/H, broadly excluding the evolution of 
oxygen; (2) the average abundances of the a-elements are consistent with those of H ii regions, 
indicating a negligible global enrichment in the disk of M33 from the epoch of the formation 
of the PN progenitors to the present time; (3) the radial oxygen gradient across the M33 disk 
has a slope of -0.031±0.013 dex kpc~^, in agreement, within the errors, with the corresponding 
gradient derived from H ii regions. Our observations do not seem to imply that the metallicity 
gradient across the M33 disk has flattened considerably with time. We report also the discovery 
of a PN with Wolf-Rayet features, PN039, belonging the class of late [WC] stars 

Subject headings: planetary nebulae: general, abundances, individual (M33 system) — galaxies: individ- 
ual(M33), evolution, abundances 

1. Introduction 



^ Observations reported here were obtained at the MMT 
Observatory, a joint faciUty of the Smithsonian Institution The galaxy M33 (NGC 598) is One of the closcst 

and the University of Arizona. Spiral galaxies of the Local Group. The measured 

distance to M33 ranges from 730 kpc (Christian 
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& Schommer 1987) to 910 kpc (Kim et al. 2002), 
with the most recent estimates of Sarajedini et 
al. (2006) and Bonanos et ah (2006) in the mid- 
range. The proximity of M33, together with its 
large angular size (optical size 53'x83', Holmberg 
1958), and its intermediate inclination {i—53°), al- 
lows detailed studies of its stellar populations and 
ionized nebulae. 

M33 is a galaxy rich in both PNc and H ii re- 
gions. The population of PNc in M33 was early 
investigated using an objective-prism survey by 
Ford (.1983) and Lequeux et al. (1987,) . The advent 
of wide-field C CD ca meras allowed deeper surveys. 
Magrini et al. 1200(f) identified 13 1 cand idate PNe. 
More recently, CiarduUo et al. (j2004l ) confirmed 
a large number of previously detected PNe and 
identified new candidates, leading to the current 
number of spectroscopically confirmed PNe: 138 
in the M33 disk and 2 in its halo. 

Due the closeness of M33 and the brightness 
of the giant H ii regions, their spectroscopy was 
obtained since the 70s. Searle (jl97ll ) presented 
spectrophotometry of eight H li regions, and fur- 
ther sp ectro scopic studies were carri ed on by 
Smith (ll975|), Kwitter & AUer (Il98l[) . and by 
Vilchez et al. (|l988f ). Several catalogs of the M33 
H II regions have been published, such as those 
by Courtes e t al. dlMZl), Calzetti et al. (Il995h . 
Wyder et al. (|l997l ). and Hodge et al. (l99^. Re- 
cently, the Local Group Census (LGC, Corradi 
& Magrini 2006) and the Local Group Survey 
(LGS, Massey et al. 2007) have provided deep 
narrow- and broad-band ^2 square degree cover- 
age of M33, disclosing a conspicuous number of 
new emission-line objects located at all galacto- 
centric distances. 

Given the wealth of information on its PN and 
H II region populations, M33 is indeed an ideal 
candidate to test chemical evolution predictions. 
The aim of the spectroscopic studies of PNe and 
H II regions is devoted both to the individual 
studies of these objects, and to draw the radial 
metallicity gradient across the disk. In particu- 
lar, PNe and H ii regions have two very differ- 
ent formation ages, and the comparison of these 
two populations provides insight on the evolution 
of the host galaxy. PN progenitors are low- and 
intermediate-mass stars (LIMS), with masses be- 
tween 1 and 8 Mq, which must have been formed 
between ~3xl0^ yr and 10 Gyr ago (Maraston 



2005). Instead H ii regions are very young. 

During their evolution, LIMS do not modify, at 
least at zcroth approximation, the composition of 
a- elements such as oxygen, neon, argon, and sul- 
fur. These elements are produced mainly from the 
nucleosynthesis of Type II supernovae and main- 
tain their capability to testify their original pres- 
ence in the interstellar cloud that gave birth to 
the PN progenitor. While there are indications 
that for extremely low metallicities both oxygen 
and neon could be modified, as it has been ob- 
served both in the SMC and in other dwarf galax- 
ies (Leisy and Dennefeld 2006; Magrini et al. 2005; 
Kniazev et al. 2008), at the metallicity of M33 
one does not expect any nucleosynthesis activity 
involving these elements (Marigo 2001). 

On the other hand, the helium, nitrogen, and 
carbon abundances measured in PNe do not cor- 
respond to those at the time of the progenitor's 
formation, since these elements are synthesized in 
LIMS. These elements hence give information on 
LIMS evolution accordingly to their initial mass 
and metallicity, and, at given metallicity, are help- 
ful to constrain the PN progenitor mass and age. 

The M33 metallicity gradient derived from a- 
element abundances has been the focus of sev- 
eral studies already, both involving PNe and H ii 
regions. The advantage of M33 with respect to 
our Galaxy, as a playground for gradients, is that 
PNe in M33 have well determined galactocentric 
distances (with relative errors within 5%) com- 
pared to the large indetermination of Galactic PN 
distances (Stanghellini et al. 2008). Magrini et 
al. (2004) using PN spectroscopy in M33 deter- 
mined the elemental abundances of 11 PNe, and 
Magrini et al. (2007b) (hereafter M07b) derived 
an oxygen radial gradient of A(0/II) / AR = - 
0.11±0.04 dex kpc~^, where R is the galactocen- 
tric distance. The sample of 11 PNe was too small 
for a definite answer to the PN gradient of M33, 
especially given the paucity of PNe observed at 
large galactocentric distances. 

The metallicity gradient of M33 using H ii 
regions has been obtained by many authors, 
with broadly differe nt results. The first stud- 
ies by Smith (Il975l). K witter & AUer (Il98l[ ). 



and Vflchez et al. (|1988I ) agreed on a steep oxy- 
gen gradient. Garnett et al. ( 1997t ) also ob- 
tained a steep gradient, with A(0/H) / AR = 
-0.11±0.02 dex kpc^^ by homogeneously com- 
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piling published data. Recent det ermin ations, 
such as th ose by Crockett et al. ( 2006[ ) , Ma - 
grini et al. (|2007a[ l. and Rosolowsky et al. (|2008h 
(hereafter M07a and RS08, respectively) seem to 
converge to a much shallower gradient, respec- 
tively deriving A(0/H) / AR = -0.012±0.011, 
-0.054±0.011, and -0.027±0.012 dex kpc- ^ The 
most recent result by Rubin et al. ( 2008[ ). based 
on Ne/H and S/H, gives -0.058±0.014 and - 
0.052±0.021 dex kpc-\ respectively. The RS08 
sample is the largest to date, and should be used 
preferentially, since the results from small samples 
might emphasize the scatter rather than the slope. 

Metallicity gradients in M33 have also been 
estimated form young giant stars (Herrero et 
al. 1994, McCarthy et al. 1995, Venn et al. 1998, 
Monteverde et al. 1997, 2000, Urbaneja et al. 2005), 
and from AGB (Cioni et al. 2008), RGB stars 
(Stephens et al. 2002, Kim et al 2002, Galletti et 
al. 2004, Tiede et al. 2004, Brooks et al. 2004, 
Barker et al. 2006), and Cepheids (Beaulieu et 
al. 2006). 

One of the most discussed questions about the 
metallicity gradient in disk galaxies is how it 
evolves with time. Chemical evolution models (see 
e.g. M07b for a review of M33 models) predict dif- 
ferent temporal behaviors of the metallicity gra- 
dient depending on assumptions such as gas in- 
flow and outflow rate, and star and cloud forma- 
tion efhciencies. Observations are needed to con- 
strain these theoretical scenarios, but so far they 
have been insufficient, especially for the old popu- 
lations. Comparing different sets of results for the 
young stellar populations, such as H ii regions, and 
the old population, such as AGB and RGB stars, 
is also delicate, since the techniques of observing 
and analyzing nebulae and stars are very different, 
each with its own collection of uncertainties. 

The idea behind the observations leading to 
this paper is to study the chemical and physical 
properties of a large number of PNe and H ii re- 
gions using the same set of observations, the same 
data reduction and analysis techniques, and iden- 
tical abundance determination methods, in order 
to avoid all biases due to the stellar vs. nebu- 
lar analysis. The chemical properties of PNe and 
H II regions will provide us with snapshots at two 
epochs in the life-time of M33, in particular of its 
metallicity gradient. In the present paper we show 
the results obtained from the PN data. We have 



obtained the first sizable sample of M33 PNe with 
uniformly derived abundances in order to build the 
first sound gradient determination from PNe. In 
a forthcoming paper we will present our own H ii 
region results, and we will compare them with the 
PN results presented here. 

This paper is organized as follows: in § [2] we 
present the observations and data reduction, in 
§ [3] we discuss the plasma diagnostics used and 
how we determine the abundances of the PNe, and 
in § |4] we describe the PN properties. In § [5] we 
compare the results with the PN populations in 
other galaxies. The chemical abundance gradients 
and their implications in the evolution of M 33 
are discussed in § [B] and in § [T] Finally in § |S] we 
summarize our results, and present the conclusions 
of this work. 

2. Observations and data reduction 

We obtained spectra of 102 PNe and 48 H ii 
regions in M33 using the MMT Hectospec fiber- 
fed spectrograph (Fabricant et al. 2005) which is 
equipped with an Atmospheric Dispersion Correc- 
tor. The spectrograph was used with a single 
setup: 270 mm~^ grating at a dispersion of 1.2 
A pixel" ^. The resulting total spectral coverage 
ranged from approximately 3600 A to 9100 A, thus 
including the basic emission-lines necessary for the 
determination of their physical and chemical prop- 
erties. The instrument deploys 300 fibers over a 
1° diameter field of view and the fiber diameter is 
~ 1.5" (6 pc using a distance of 840 kpc to M 33). 
The 102 PNe w ere se lected from the catalog of 
CiarduUo et al. ( 2004h . including four new PNe 
discovered therein at large galactocentric radii. 

In Table [T] we present the list of the observed 
PNe. Column (1) give s the identification number 
from Ciardullo et al. (|2004f ) except for PNe 153 
to 156 which are from LGC observations (Corradi 
& Magrini 2006); columns (2) and (3) give the 
equatorial coordinates, RA and DEC at J2000.0, 
and column (4) gives the [O ill] magnitude from 
Ciardull o et a l. (|2004[ ) following the definition by 
Jacoby (Il989l ). 

The H II regions of our sample were chosen ei- 
ther based on their [O ill] brightness, or among 
those whose chemical abundances had already 
been measured by others so we can use them as 
a control sample. In this paper we only use the 
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Table 1 
Observed PNe 



Id 


RA J2000.0 DEC 


M [OIII] 


(1) 


(2) 


(3) 


(4) 


PNOOl 


1:32:09.04 


30:22:05.700 


23.51 


PN002 


1:32:26.54 


30:25:49.800 


23.53 


PN003 


1:32:38.03 


30:24:00.603 


21.82 



Note. — (1) Identification number from Ciardullo 
et al. II2004I) . except for PNo 153 to 156, which 
are from LGC observations (Corradi & Magrini 
2006); (2), (3) galactic coordinates at J2000.0; (4) 
[O III] mag nitude following the definition by Ja- 
coby (119891) . Table 1 is published in its entirety 
in the electronic edition of the Astrophysical Jour- 
nal. A portion is shown here for guidance regarding 
its form and content. 



data of the H ii regions to test our procedures by 
comparing our fluxes with those in the hterature. 

We used a large number of fibers to take sky 
spectra for sky-subtraction. In order to opti- 
mize the sky-subtraction we selected low diffuse- 
emission areas on the face of M33 from the INT 
Ha and [O iii]-continuum frames (Magrini et al. 
2000). Since the continuum from the PN central 
star can not be detected at the M33 distance any 
continuum emission in the spectra has to be aris- 
ing from an underlying unresolved stellar popu- 
lation. However, the continuum emission in our 
spectra is in most cases negligible, and therefore 
the Balmer absorption lines coming from the stel- 
lar background are unimportant. 

Five 1800 s exposures were taken on the night of 
October 13, 2007, and three additional 1800 s ex- 
posures of the same field were obtained on Novem- 
ber 12, 2007. The airmass during the observa- 
tions ranged between 1.07 and 1.4 during the night 
of October 13, and from 1.3 to 1.5 on November 
12. The seeing was 1.2 " and 1.7", respectively. 
Several dome-flat and sky-flat exposures were ob- 
tained during the nights of observations to perform 
the data reduction. Arc exposures with the cali- 
bration lamp He-Ne-Ar were taken for wavelength 
calibration. 

The spectra were reduced using the Hectospec 
package. All observations were bias subtracted, 
overscan corrected, and trimmed. The science ex- 
posures were flat-fielded and combined together 
to eliminate cosmic rays, and the one-dimensional 
spectra were extracted and wavelength calibrated. 



The relative fiux calibration was done observing 
the standard star Hiltm600 (Massey et al. 1988) 
during the nights of October 15 and November 
27. The standard star was observed with airmass 
^1.1. The emission-line fluxes were measured with 
the package SPLOT of IRAjH. Errors in the fluxes 
were calculated taking into account the statisti- 
cal error in the measurement of the fluxes, as well 
as systematic errors of the flux calibrations, back- 
ground determination, and sky subtraction. 

The observed line fluxes were corrected for the 
effect of the interstellar extinction using the ex- 
tinction law of Mathis (1990) with i?y=3.1. 

We derived c(H/3), the logarithmic nebular ex- 
tinction, by using the weighted average of the 
observed-to-theoretical Balmer ratios of Ha, H7, 
and H<5 to H/3 (Osterbrock & Ferland 2006). 

Table 2 gives the results of our line measure- 
ments and extinction corrections. Column (1) 
gives the PN name; column (2) gives the nebular 
extinction coefficient c(H/3) with its error; columns 
(3) and (4) indicate the emitting ion and the rest- 
frame wavelength in A; columns (5), (6), and (7) 
give the measured line fiuxes (F^), their absolute 
errors (A(Fa), and finally the extinction corrected 
fiuxes (Ia). Both Fa and I a are normalized to 
H/3=100. 

In order to confirm the goodness of our spec- 



^IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Uni- 
versities for Research in Astronomy (AURA) under coop- 
erative agreement with the National Science Foundation 
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Table 2 

Observed and de- reddened fluxes. 



Id 


c(H/3) 


Ion 


A (A) 


Fa 


A(Fa) 


Ia 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




U.OZO±U.UzD 


HI 


4340 


47.6 


3.7 


52.4 






HcII 


4686 


58.0 


4.4 


60.1 






HI 


4861 


100.0 


4.8 


100.0 






[OIII] 


4959 


304.2 


7.4 


298.9 






[OIII] 


5007 


917.8 


12.2 


894.0 






HI 


6563 


362.2 


8.0 


289.0 


PN002 


0.343±0.030 


HI 


4340 


42.3 


3.9 


46.9 






HI 


4861 


100.0 


4.8 


100.0 






[OIII] 


4959 


304.1 


7.5 


298.5 






[OIII] 


5007 


953.0 


12. 


927.04 






HI 


6563 


368.3 


8. 


289.70 


PN003 


0.460±0.014 


loiil 


3727 


101.8 


1.6 


137.2 






HI 


3835 


7.7 


0.8 


10.1 






Hel 


3889 


9.1 


0.8 


11.7 






[NeIII]/HI 


3968 


7.6 


0.6 


9.6 






HI 


4100 


17.3 


0.8 


21.1 






HI 


4340 


37.5 


1.2 


43.1 






[OIII] 


4363 


7.4 


0.7 


8.5 






Hoi 


4471 


4.4 


0.7 


4.9 






HI 


4861 


100.0 


1.8 


100.0 






[OIII] 


4959 


178.5 


2.0 


174.2 






[OIII] 


5007 


540.1 


3.0 


520.2 






[Nil] 


5755 


3.7 


0.7 


3.0 






Hoi 


5876 


22.4 


1.0 


18.0 






[Nil] 


6548 


43.5 


1.0 


31.6 






HI 


6563 


384.8 


2.4 


280.0 






[Nil] 


6584 


128.2 


1.6 


92.7 






Hoi 


6678 


5.2 


0.7 


3.7 






[SII] 


6717 


7.2 


0.6 


5.1 






[SII] 


6731 


10.3 


0.8 


7.3 






Hoi 


7065 


9.7 


0.8 


6.5 






[Arlll] 


7135 


14.8 


0.8 


9.9 



Note. — (1) PN name; (2) nebular extinction coefficient c(H/3) with its error; 
(3) emitting ion; (4) rest-frame wavelength in A; (5) measured line fluxes; (6) 
absolute errors on the measured line fluxes; (7) extinction corrected line fluxes. 
Both F;^ (5) and (7) are expressed on a scale where H/3— 100. Table [2] is 
published in its entirety in the electronic edition of the Astrophysical Journal. 
A portion is shown here for guidance regarding its form and content. 
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troscopic calibration we compare our emission-line 
flux measurements of both H ii regions and PNe 
with previously published fluxes of the same ob- 
jects. In Figure 1 we plot our measured fluxes 
against the ones from the literature, where the 
two dashed lines mark differences of ±0.15 dex 
between the two sets. We found a good agree- 
ment between the sets, especially for bright emis- 
sion lines, giving us confidence of a sound spectral 
calibration. 

In Figure [2] we show, as an example of the data 
quality, the spectra of three PNe with different 
excitation. 

3. Plasma Diagnostics and Abundances 
3.1. Electron Densities and Temperatures 

We have used the extinction-corrected intensi- 
ties to obtain the PN electron densities and tem- 
peratures. In order to calculate the electronic 
densities we used the doublet of the sulfur lines 
[S li]AA6716,6731, while for the electron temper- 
atures we used the ratios [O iii]A4363/(A5007 -|- 
A4959) and [N ii]A5755/(A6548 + A6584). We 
performed plasma diagnostics by using the 5-level 
atom model included in the nebular analysis pack- 
age in IRAF/STSDAS (Shaw & Dufour 1994). 

We have also used the standard forbidden line 
diagnostics IRAF routines in nebular to deter- 
mine the electron temperatures. We have deter- 
mined the low- and medium-excitation tempera- 
tures from the [O III] and the [N ll] line ratios, 
respectively (see also Osterbrock & Ferland 2006, 
§5.2). The temperature uncertainties have been 
estimated by formal error propagation of the ab- 
solute errors on the line fluxes (see also Table 4). 
The average relative uncertainties in the determi- 
nation of Te[0 III] and Tc[N II] are of the order 
of 5% and 15%, respectively. 

The [O III] and [N ll] temperature diagnostics 
are available, respectively, for 32 and 8 PNe of our 
sample. Only for six PNe we can use both temper- 
ature diagnostics. In order to improve the size of 
the PN sample with temperature determinations 
we looked for correlations between the electron 
temperature and various diagnostics of the neb- 
ular excit ation, similarly to the approach followed 
by Kaler ( 1986f ) for a Galactic PN sample. 

In Figure [3] we plot the [O ill] electron tem- 



o 



6£1 
O 




1 1 

log(F,/F„^), ref 



Fig. 1. — Emission- line flux measurements of both 
H II regions (filled circles) and PNe (empty cir- 
cles) are compared to previously published fluxes. 
The fluxes of the H ii reg ions have been taken 
from Kwitter & AUer (|l98l|); Vilchez et al. (198^; 
Crockett et al. (|2006l ). M07a ; RS0 8, and the PN 
ones are from Magrini et al. (j2003l ). 



2. 5x10" 




20 40 60 80 

I (A4686) 



Fig. 3. — [O III] electron temperature vs. Ia4686 
line intensity, scaled for Ih/3=100 and corrected for 
extinction. The solid line is the fit done with the 
fitexy routine, shown in Eq. (1). 
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2xlO-if 




4000 5000 6000 



wavelength (A) 

Fig. 2. — Spectra of three PNe with different excitation selected from our sample in order to show the data 
quality. 
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perature against the intensity of the Ia4686 line, 
scaled for Ih/3=100 and corrected for extinction. 
Both quantities are plotted with their formal er- 
ror bars. Planetary nebulae in this plot are those 
hot enough for helium to be doubly ionized. We 
find a clear correlation between the [O ill] electron 
temperature and the intensity of the Ia4686 Ihie, 
an effect that is expected as the central stars heat 
the medium-high excitation nebulae. In order to 
use this correlation in our temperature calculation 
we have applied the routine fitexy in Numerical 
Recipes (Press et al. 1992) to fit the relation be- 
tween the two variables, taking into account their 
errors, and minimizing x^- 
For Ia4686 >0 we find: 

Tc [O III] - (63 ± 8) Ia4686 + (11, 000 ± 200), (1) 

where To[0 III] is given in K and Ia4686 is scaled 
to 11/3=100. Note that using a different fitting 
routine, such as sixlin (Isobe et al. 1990), we de- 
termine temperatures that agree with the ones of 
Eq. (1) at the 1% level. We infer that we are able 
to estimate the temperatures to better than 10% 
with the above formulation. 

If Ia4686=0 we assum e a constant value of To[0 
III], as in Kaler (il986} ). The average To[0 III] for 
Ia4686=0 is 12100 ±1500 K, which is compatible 
with the continuity of Eq. (1). Note that both this 
value and Eq. (1) give temperatures that ar e con - 
sistently higher than those found by Kaler ( 1986f ) 
in Galactic PNe, an effect that is probably due to 
the lower metallicity that characterizes th e M33 
stellar populations (see Stanghellini et al. ( 20031 ) 
for a discussion on how metallicity affects the neb- 
ular temperatures). 

In Figure [4] we show Te [N II] against Ia4686 , 
where both quantities are plotted with their for- 
mal errors. The line represents the mean value of 
Te[N II]= 13600 ±2100 K, which has been derived 
using the (inverse squared) relative temperature 
errors as weights. Once again, while the physi- 
cal behavior of the te mperatures is similar to that 
found in Kaler (jl986|), the actual values are higher, 
possibly reflecting the lower efficiency of the cool- 
ing process at low metallicity. For Ia4686 >0 we 
adopt Te[N II] = 13,600 K. 

If Ia4686=0, we expect To[N II] to decrease with 
the stellar temperature. The only way to frame 
this decline is to use the ratio between the fluxes 
of the singly-ionized to the doubly-ionized oxygen. 



By correlating the oxygen line strengths, corrected 
for extinction, to the electron temperatures for the 
PNe with measured [N ii] lines we find: 



To [N II] = (ll,700±2800)-(2800±4000)log ^ AA3727 29 ^ 

Ia4959 ± Ia5007 

(2) 

Note that in our spectra the AA3727 — 29 
lines are always blended. If the oxygen lines 
are not available, and Ia4686==0 we assume Tc[N 
II] = 14, 100 ±2800 K, which, in this case, is the 
weighted mean of the temperatures. 

3.2. Chemical abundances 

We computed the ionic abundances using the 
TieftiiZar analysis package in IRAF/STSDAS (Shaw 
& Dufour 1994). The elemental abundances were 
then determined by applying the ionization correc- 
tion factors (ICFs) following the prescriptions by 
Kingsburgh & Barlow (1994i) for the case where 
only optical lines are available. As discussed in 
the previous §, the [O ill] A4363 emission line was 
measured with a sufficiently high signal to noise 
ratio in 32 PNe, affording the direct determina- 
tion of the Tc[0 III]. We could also directly mea- 
sure To[N II] in 8 PNe. For the other targets we 
adopted the temperatures derived from the rela- 
tions we have described in the previous Section. In 
the abundance analysis we have used To[N II] for 
the calculation of the N+, 0+, S+ abundances, 
while Te[0 III] was used for the abundances of 
02+, S2+, Ar2+, He+, and He2+. 

We calculated the abundances of He i and H e 11 
using the equations of Benjamin et al. ( 19991 ) in 
two density regimes, i.e. Ue >1000 cm~'^ and 
<1000 cm^'^. The Clegg's coUisional populations 
were taken into account (Clegg 1987). Due to 
the effect of double collisions not being properly 
corrected, the He+/H ionic abundance from the 
A7065 emission-line is quite different from the 
ionic abundance computed from AA4471, 5876, 
and 6678. Wc thus omit the abundance derived 
from the former line from the average. The com- 
puted individual ionic and elemental abundances 
of each PN are shown in Table 3, where column (1) 
gives the PN identification name; columns (2) and 
(3) show the plasma diagnostic and abundances 
available for each PN, where in column (2) we la- 
bel each diagnostic, and in column (3) we give the 
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Table 3 

Plasma Dl^gxostics and Abundances 



Id 






(1) 


(2) 


(3) 



PNOOl 



Tf^fTHTT^t 
J. t3 1 11 1 I 


14650^ 


±C 1 iN 11 ] 


13600^ 


HcII/H 


063 


Hc/H 


0.063 


OIII/H 


7.586 10~^^ 


ICF(O) 


1 000 


O/H 


7.586 10^"^ 


PN002 




To(OIII) 


14650*' 


To(NII) 


13600*' 


OIII/H 


7.884 10"*"^ 


ICF(O) 


1.000 


O/H 


7.884 10-°^ 


PN003 




Te(OIII) 


14000" 


Te(NII) 


15200'' 


Te(OIII) 


12060*' 


To(NII) 


14120*' 


Nc(SII) 


2540 


Hcl/H 


0.125 


Hc/H 


0.125 


OII/H 


1.590 10""^ 


OIII/H 


6.500 10"°^ 


ICF(O) 


1.000 


O/H 


8.090 10-°^ 


NII/H 


7.168 10"°" 


ICF(N) 


5.087 


N/H 


3.646 10-°^ 


Arlll/H 


4.47 10-°'' 


ICF(Ar) 


1.245 


Ar/H 


8.359 10-°'' 


SII/H 


1.736 10-°'' 


ICF(S) 


1.276 


S/11 





^computed from electron tempera- 
ture diagnostic lines. 

"^derived from the relations in §3. 



Note. — (1) identification name: (2) 
laljcl of each plasma diagnostic and 
abundances available; (3) relative 
value obtained from our analysis. Ta- 
ble 3 is published in its entirety in the 
electronic edition of the Astrophysical 
Journal. A portion is shown here for 
guidance regarding its form and con- 
tent. 
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Table 4 
Typical Errors in dex 



Magnitude 


Tc 


A(Ho/H) 


A(0/H) 


ACN/H') 


A(Nc/H) 


A(Ar/H) 


ACS/H) 


N PNc 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


20.60-21.00 


1 


0.03 


0.025 


0.06 


0.05 


0.06 


0.09 


8 




c 




0.17 


0.40 


0.43 


0.34 


0.43 


1 


21.01-21.50 


1 


0.04 


0.03 


0.06 


0.06 


0.11 


0.11 


9 




c 


0.09 


0.20 


0.34 


0.43 


0.30 


0.30 


8 


21.51-22.00 


1 


0.04 


0.03 


0.04 


0.08 


0.09 


0.09 


5 




c 


0.08 


0.21 


0.34 


0.43 


0.34 


0.30 


8 


22.01-22.50 


1 


0.08 


0.04 


0.07 


0.15 


0.17 


0.11 


6 




c 


0.11 


0.21 


0.40 


0.43 


0.40 


0.30 


16 


22.51-23.00 


1 


0.13 


0.09 


0.16 


0.16 


0.17 


0.11 


4 




c 


0.17 


0.20 


0.40 


0.43 


0.43 


0.30 


9 


23.51-24.00 


1 


0.09 


0.07 


0.09 


0.13 


0.18 


0.16 


1 




c 


0.21 


0.23 


0.40 


0.43 


0.43 


0.30 


12 


24.01-24.50 


c 


0.11 


0.21 


0.34 


0.43 


0.43 


0.43 


5 


24.51-25.00 


c 




0.45 










2 


25.01-25.50 


c 




0.5 










1 



Note. — (1) Magnitude range; (2) method of the electron temperature determination, where I means that 
the temperature has been measured based on an emission line and c means that we used the correlations 
described in §3; (3-8) typical errors in dex of the total abundances; (9) number of PNc available in each bin. 



relative value obtained from our analysis. 

The formal errors on the ionic and total abun- 
dances were computed taking into account the un- 
certainties in the observed fluxes, in the electron 
temperatures and densities, and in c(H/3). The 
errors (in dex) of the final abundances are given 
in Table ID Typical errors are given for PNe in 
different [O ill] magnitude ranges, taking also into 
account whether the electron temperatures were 
derived from emission-line diagnostic or from the 
relations described in § 13.11 In the latter case, we 
assumed a percentage error on the electron tem- 
perature of ~10%, which is the dispersion of the 
relation used to derive the temperature. The first 
column of Table S] gives the magnitude range; col- 
umn (2) gives the method used for the electron 
temperature determination, where I means that 
the temperature has been measured based on an 
emission line and c means that we used the cor- 
relations described in §3; columns (3) to (8) give 
the typical errors in dex of the total abundances. 
The last column gives the number of PNe on each 
bin. 

4. The PN population 

Most of the PNe discovered in M33 belong to 
its disk, with only two PNe having radial velocities 
compatible with the halo population (CiarduUo et 



2x10" 



^^1.5x10" 



10'' 



20 40 60 80 



Fig. 4.— Same as Figure 3 but for the Te[N II] 
electron temperature. The solid line represents the 
mean value of To[N II]= 13600 ±2100 K, which has 
been derived using the relative temperature errors 
as weights. 



10 



Table 5 

Average chemical abundances and abundance ratios 



Sample 


Ho/H 


O/H (10-'') 


N/H (10-") 


Ne/H (IQ-'') 


Ar/H (10-**) 


S/H (IQ-i^) 


N/O 


Ne/O 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


M33 PNo, Disk (91) 


0.118±0.074 


1.96±1.27 


1.09±1.54 


4.09±3.08 


1.19±0.57 


3.34±1.71 


0.40±0.32 


0.17±0.06 


M33 PNc, non-Type I (72) 


0.114±0.083 


1.73±0.93 


0.39±0.19 


3.65±2.11 


1.12±0.50 


4.46±2.30 


0.18±0.06 


0.17±0.06 


M33 PNc, Type I (19) 


0.135±0.027 


2.70±1.88 


1.66±1.88 


4.83±4.25 


1.38±0.68 


6.64±3.80 


0.58±0.33 


0.16±0.04 


M33 PNo, Halo (2) 


0.120: 


2.62±0.31 


0.84: 


3.47±0.82 


1.24±0.28 


6.45: 


0.31: 


0.13±0.03 


M33 HII regions 


O.lOliO. 015(a) 


2.04±0.75(b) 


1.28±0. 57(a) 


4.24±2.61(c) 


1.31 ±0. 45(a) 


5.85±2.28(a) 


0.06±0. 02(a) 


0.20±0. 06(e) 


Solar value(d) 


0.085±0.02 


4.57±0.04 


0.60±0.09 


6.9±1.0 


1.51±0.30 


13.8±2.0 


0.13±0.10 


0.15±0.05 


LMC PNe 


0.103±0.026 


2.32±1.65 


1.48±1.75 


4.04±3.60 


1.14±0.72 


3.46±8.88 


0.87±1.15 


0.17 ±0.09 


SMC PNe 


0.113±0.022 


1.05±0.46 


0.28±0.33 


1.77±1.32 


0.59±0.59 


4.80±6.57 


0.28±0.50 


0.17 ±0.08 


Galactic PNe 


0.123±0.042 


3.53±1.95 


2.44±3.46 


9.68±7.98 


1.26±1.24 




0.67±0.82 


0.25±0.10 



Note. — (1) the PN sample; (2-7) average eliemical abuiidanees by number with their rms uncertainties; (8-9) <N/0> and <Ne/0>, the average of N/O and Ne/O 
values for each object. 

^Chemical abundances computed using the sample of M07a, extended to all literature data with direct electron temperature measurement. 

^O/H is calculated with the cumulative sample of M07a, which includes also all previous O/H with direct electron temperature measurement, and RS08. 

"^Ne/H and Ne/O are from Crockett et al. ( 120061) . 

*^Solar value from Asplund ct al. (|2005f) . 

'Value from a single object. 



al. 2004). Therefore, we have the opportunity to 
analyze in detail a pure disk PN population by ex- 
cluding from our analysis the two suspected halo 
PNe (namely PN067 and PN024). We have also 
excluded from our analysis the PN039 since its 
physical conditions do not allow us to obtain abun- 
dance diagnostics. PN039 has a WR- nucleus and 
will be discussed in detail in the Appendix. 

4.1. The He/H vs N/O diagram 

The plot of the N/O ratio vs the He/H one is a 
classical diagnostic diagram used to discriminate 
PNc of different types, i.e. with different mass 
progenitors. The N/O ratio provides information 
about the stellar nucleosynthesis during the AGB 
phase of LIMS. On the one hand, nitrogen is pro- 
duced in AGB stars in two ways: by neutron cap- 
ture, during the CNO cycle, and by hot-bottom 
burning. Hot-bottom burning produces primarily 
nitrogen but occurs only if the base of the con- 
vective envelope of the AGB stars is hot enough 
to favor the conversion of ^^C into ^"^N. Thus ni- 
trogen is expected to be mostly enriched in those 
PNe with the most massive progenitors, i.e. with 
turnoff mass larger than '^3 Mq (van der Hoek 
& Groenewegen 1997; Marigo 2001). The He/H 
ratio gives also an indication of the initial mass 
of the progenitor, the nebula is enriched progres- 
sively for more massive stars, it reaches a plateau 
between 3 and 4 M0, and then it increases again 
toward the higher masses (Marigo 2001). 

In order to make a selection of the PNe with 
high-mass progenitors, in Figure [5l we plot the 
M33 disk PNe on the H e-N/0 plane. Peimbert & 
Torres-Peimbert (|l983l ) found that Galactic PNe 
that were nitrogen- and helium-enriched also had 
bipolar shape, and were located closer to the 
Galactic plane. It looked like these PNe formed 
a younger population, which they defined as those 
with log (N/O) > -0.3 and Hc/H>-0.125, and 
calle d the m Type I PNe. Kingsburgh & Bar- 
low (|l994l ) analyzed a much larger sample and re- 
define the Galactic Type I PNe as those having 
log(N/O)>-0.1, with 18% of the P Ne in their sam- 
ple being of Type I. Dopita (|l99ll ) analyzed LMC 
PNe and noted that the Type I definition needed 
a revision to allow for the lower metallicity of the 
LMC with respect to the Galaxy, setting the limit 
of Type I PNe to log(N /O)>-0.5 for the LMC. 
Leisy & Dennefeld ( 2006[ ) basically confirmed this 



limit. Magrini et al. (|200J) also noted that the 
definition of Type I PNe depends on the metallic- 
ity, because the amount of nitrogen that can be 
produced by hot bottom burning is dependent on 
the amount of carbon present, and also because 
the oxygen abundance depends on the metallicity 
of the galaxy. Since the oxygen metallicity of M33 
is very similar to that of the LMC, we use the Type 
I limits as in Dopita (1993), namely. Type I PNe 
are those with log(N/O)>-0.5, independent of he- 
lium abundance. With this definition we found 19 
PNe being Type I in the disk of M33 (-20% of the 
whole disk population) and one in the halo. 

4.2. The Ne/H and S/H vs O/H diagrams 

The study of the chemical evolution of galaxies 
needs strong constraints regarding the past com- 
position of the ISM, and PNe can supply them, 
in particular with their soundly-determined oxy- 
gen abundances. The questions is whether oxygen 
(and neon) is modified in LIMS while in their AGB 
phase. 

A possible way to verify this is to study the 
relationship between neon and oxygen. These el- 
ements derive both from primary nucleosynthesis, 
mostly in stars with M>10 M©. If the O/H and 
Ne/H abundances are really independent of the 
evolution of the PN progenitors through the AGB 
phase, a tight correlation between their abundance 
should be observed. 

In Figure [HI we show O/H against Ne/H for the 
55 M33 PNc where both abundances were avail- 
able. The slope of the correlation is close to unity, 
0.90±0.11, with a correlation coefficient Rxy=0.81, 
pointing at a locked variation of these t wo ele- 
ments. Recent results from Wang & Liu (2008|) 
indicate that oxygen and neon could be manufac- 
tured with similar yields also in LIMS, but only at 
very low metallicities, 12 -I- log(0/H) < 8. These 
results are however still not explained by cur- 
rent nucleosynthesis theories that predict different 
channels, and thus different yields, for Ne and O 
production in LIMS at low metallicity (cf., e.g., 
Karakas & Lattanzio 2003, Marigo et al. 2003). 
Our findings confirm that in M33, as in LMC, and 
the Galaxy, there is no evidence of enhancement 
of oxygen and neon in PNe. 

Similar correlations are expected also for sul- 
phur versus oxygen, since sulphur is manufactured 
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by massive stars. In Figure [7] we plot the sul- 
fur versus oxygen abundances of the 38 M33 PNe 
where these abundances are available. Their rela- 
tionship has a slope of 0.97±0.22 and a correlation 
coefficient of 0.51, showing moderately good corre- 
lation. We refrain from showing and analyzing the 
argon to ox ygen relation. As seen in Stanghellini 
et al. (|2006l ). these quantities might not correlate 
due to their different origin. Even if both ele- 
ments are manufactured in massive stars, argon is 
synthesized in very different a -element processes 
than oxygen, thus a lockstep of these elements is 
not completely expected. 

5. Comparison of abundances of different 
stellar populations, and of PN popula- 
tions of other galaxies 

The average chemical abundances of our M33 
PNe are reported in Table [5] where the first col- 
umn describes the selected sample, columns (2) 
to (7) give the average total abundance of He/H, 
0/H, Ne/H, Ar/H, and S/H, whereas the last 
two columns, (8) and (9), give the mean values 
of N/0 and Ne/0. The average abundances were 
computed excluding those abundances from up- 
per limit detections. When the emission lines 
for temperature diagnostics were not detected, we 
adopted the electron temperatures discussed in 
§3.11 Their larger uncertainties were taken into 
account by larger formal errors in the total chem- 
ical abundances as described in Tabled) 

We group the PNe in M33 into four popula- 
tions: i) the whole PN disk population; it) the 
disk population excluding those PNe which are be- 
lieved to be Type I; the Type I PNe, and iv) 
the ha lo PN population as classified by CiarduUo 
et al. (|2004j ). For the sake of comparison we add 
to the Table the average chemical abundances of 
the M33 H ii regions, the solar abundances, and 
the average abundances of the PNe in LMC, SMC 
(Stanghellini 2008) and the Galaxy (Stanghellini 
et al. 2006). Note, however, that the population of 
the Type I and non-Type I classes remains some- 
what uncertain, since about a dozen PNe could 
belong to either class based on their formal abun- 
dance errors. However, the chemical abundances 
and distributions of Type I and non-Type I PNe 
in M33 are very similar, so that the displacement 
of a few PNe from a group to the other would not 



0.5 



-0.5 



o 



■1.5 




10.8 11 11.2 

12 + log(He/H) 



11.4 



Fig. 5.— log(N/0) versus 12-hlog(He/H). Type 
I PNe and non-type I PNe are plotted as filled 
circles and empty circles r espect ively, according to 
the definition of Dopita ( 199ll ). The continuous 
line mark log(N/O)>-0.5 that, independently of 
helium abundance, defines the area of Type I PNe. 




8.5 9 
12 + log(0/H) 



Fig. 6. — The relationship between oxygen and 
neon abundances. Symbols are as in FigureO The 
continuous line is the weighted least square fit to 
the complete sample of PNe. 
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change our results significantly. 

An inspection of Table 5 affords several inter- 
esting clues on the chemical evolution of galaxies. 
First, if we compare the average abundances of 
the a-elements of the M33 disk PNe with those of 
the H II regions, we do not note significant enrich- 
ment in the H ii regions, leading to the conclusion 
that these elements are not substantially modified 
in the LIMS evolution. We substantially found 
a flat age-metallicity relationship in M33 all the 
way to the time of PN formation. Second, we in- 
fer that the metallicity of the PNe in M33, based 
on oxygen abundances, is sub-solar. If we consider 
the entire sample of disk PNe in M33 we see that 
the averages of a-elements are virtually identical 
to those of PNe in the LMC, with PN metallicity 
about 1/2 that of the Galaxy, and twice the SMC 
on average, both in oxygen and neon. 

From Table 5 we see that the average Ne/0 ra- 
tio in the general population of M33 PNe is iden- 
tical to that of the LMC, but lower than the mean 
Ne/0 in Galactic PNe. This finding is right on th e 
mark with the recent results of Wang & Liu (2008) . 
In fact Wang & Liu (fiooal found that the Ne/0 
ratio increases with increasing oxygen metallicity 
both in PNe and in H ii regions. This suggests 
a different enrichment history of neon and oxygen 
in the ISM and thus probably different produc- 
tion mechanisms of these two a-elements in mas- 
sive stars, in agreement with c urrent theoretical 
calculations by Kobayashi et al. While we 

should expect the solar Ne/0 ratio to be consistent 
with that of Galactic disk PNe and H ii regions, 
we note that the Asplund et a l.'s (ioosh value, 
as well as Grevesse & Sauval 's (1998^ is slightly 
discrepant in this scenario. Both PNe and H ii re- 
gions in the Galactic disk give a consistent Ne/0 
ratio of ~0.25, higher than the solar value. Our 
results on M33 PNe seem thus to be consistent 
with the recent suggestions of a revision of the so- 
lar Ne/0 ratio and the absolute neon (e.g. Wang 
& Liu 2008, Rubin et al. 2008). 

The N/0 average shown in Table 5 for M33 
PNe is lower than the Galactic and higher than the 
SMC values. The comparison to the LMC average 
has limited importance, giving the large range of 
LMC N/0 values. 



6. The abundance gradients 

In Figure [H we show the oxygen abundance 
as a function of galactoccntric distance for our 
sample of disk PNe. The galactoccntric distances 
were computed by adopting a distance to M33 



840 kpc, as determined by Frecdman et al. (jl991r ). 
an average and well-accepted distance estimate. 
A weighted linear least-square fit to the complete 
disk sample gives a gradient: 



12-|-log(0/H) = -0.031(±0.013)Rgc+8.44(±0.06), 

(3) 

where Rgc is the de-projected galactoccntric 
distance in kpc, computed assuming an inclination 
of 53°, and a position angle of 22°. In Figures [U 
and [TUl we plot the radial metallicity gradients of 
Ne/H and S/H. The solid lines are the weighted 
linear fits to the complete sample of disk PNe (first 
row of Table [6] for each element) . In all these fig- 
ures the two symbols refer to Type I (filled circles) 
or non-Type I (empty circles) PNe. 

In Table [6] we show the slopes (column 2) and 
zero-points (chemical abundances in the centre of 
the galaxy, column 3) of the metallicity gradients 
considering different samples of PNe. The sample 
used and the number of PNe included are given 
in columns (4) and (5) respectively. For each el- 
ement, the first row report the gradient obtained 
considering the whole sample of disk PNe (thus 
excluding the two possible halo PNe), the second 
row gives the gradient computed using non-Type I 
PNe, and the third row gives the gradient obtained 
for Type I PNe. 

The gradients are computed with a weighted 
least mean square fit. For each element the slope 
of the three different samples are consistent within 
the errors. Note that the metallicity gradient of 
an old stellar population could be affected by the 
radial migration of stars during cosmic times. The 
present-time location of a PN could not necessar- 
ily be the place where it was born. Radial mixing 
of stars is believed to be due to several mecha- 
nisms, among them the diffusion of stars on their 
orbits because of various irregularities in the galac- 
tic potential (cf. Wielen et al. 1996), to the pas- 
sage of spiral patterns (cf. Minchev & Quillcn 
2006), to changes in the angular momentum due 
to non-asymmetric forces due to molecular clouds 
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(Spitzer & Schwarzzchild 1953). Upper limits to 
the radial migration rate were estimated by De Si- 
mone et al. |20O4) and by Haywood (2008), giving 
similar values around 1.5-3.7 kpc Gyr^^, limited 
to a radial region ±2 kpc fro m the b irth place. The 
model of Sellwood & Binney (|2002l) show that spi- 
ral waves in galaxy disks churn both the interstel- 
lar medium and the stars, affecting their metallic- 
ity gradients. The effect of stirring of the entire 
disk due to the spiral waves is a flattening of the 
metallicity gradient. 

However, the fact that metallicity gradients sur- 
vive and are well observed in disk galaxies seem 
to indicate that the effectiveness of the migration 
processes is moderate. 

7. Discussion 
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The metallicity gradient of PNe allows to an- 
alyze how was the metal distribution in the past 
epochs of M33. In M33 we have to possibility to 
compare the metallicity gradient of PNe with a 
good number of other measurements of abundance 
gradients, in particular from H ii regions. The ad- 
vantage of comparing PNe to H ii regions is the 
similarity of their emission-line spectra in spite of 
their different evolutionary states that allows to 
use the same observation techniques, analysis, and 
abundance determinations. 

In order to compare our PN results with the 
leading samples of M33 H ii regions we build the 
largest, most reliable sample available to date. As 
already noted by e.g. RS08 and M07a, the size 
and quality of the sample is fundamental to de- 
rive reliable metallicity gradients. In H ii regions 
there is a substantial intrinsic scatter of 0.11 dex 
in the metallicity at any given distance from the 
M33 center, which imposes a fundamental limit on 
the accuracy of gradient measurements that rely 
on small samples of objects (RS08). Also, in the 
past, only giant H II regions have been studied, 
and those might have a steeper metallicity gradi- 
ent than small and compact H ii regions. 

We have re-calculated the H ii regions metallic- 
ity gradients considering a cumulative sample, in- 
cluding that by RS08, consisting of 61 H ii regions 
with accurate abundance determinations, and that 
by M07a that included new and literature spec- 
troscopy performed to a set of 28 H ii regions. 
The cumulative sample of 89 H ii regions gives. 



Fig. 7. — The relationship between oxygen and 
sulphur abundances. Symbols are as in Figure [5] 
The continuous line is the weighted least square 
fit to the complete sample of PNe. 




Fig. 9. — The radial gradient of neon abundance. 
Symbols are as in Fig.[5l The continuous line is the 
weighted least square fit to the complete sample 
of disk PNe. Slopes and zero-points are shown in 
Table [3 
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Fig. 8. — The radial gradient of oxygen abundance. Symbols are as in Fig. [S] The continuous line is the 
weighted least square fit to the complete sample of disk PNe. Slopes and zero-points are shown in Table [6] 



16 



Table 6 

Chemical abundance gradients 



Element 


Slope 


Zero 


sample 


# PNc 


(1) 


(2) 


(3) 


(4) 


(5) 


O/H 


-0.031±0.013 


8.44±0.06 


All PNe 


91 


O/H 


-0.030±0.013 


8.41±0.06 


non-type I PNe 


72 


O/H 


-0.039±0.033 


8.56±0.15 


type I PNc 


19 


No/H 


-0.037±0.018 


7.75±0.07 


All PNe 


55 


No/H 


-0.051±0.019 


7.76±0.08 


non-type I PNe 


37 


No/H 


-0.025±0.037 


7.76±0.17 


type I PNo 


18 


S/H 


-0.033±0.019 


6.77±0.08 


All PNe 


38 


S/H 


-0.027±0.020 


6.67±0.09 


non-type I PNe 


19 


S/H 


-0.041±0.031 


6.87±0.14 


type I PNc 


19 



Note. — (1) The chemical element for which the gradient is computed; 
(2) slope of the abundance gradient; (3) zero-point (chemical abundances 
in the centre of the galaxy); (4) sample of PNe; (5) number of PNe 
in each sample. For each element, the first row report the gradient 
obtained considering the whole sample of disk PNe (thus excluding the 
two possible halo PNe), the second row gives the gradient computed 
using non-Type I PNe, and the third row gives the gradient obtained of 
Type I PNe. 




Fig. 10. — The radial gradient of sulphur abun- 
dance. Symbols are as in Fig. [5l The continuous 
line is the weighted least square fit to the com- 
plete sample of disk PNe. Slopes and zero-points 
are shown in Table [G] 



Fig. 11. — The slopes of the radial elemental gra- 
dients in M33. The continuous line is the weighted 
average slope, and the two dashed line give the rms 
a of this value. 
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12+log(0/H) = -0.032(±0.009)xRgc+8.41(±0.04). 

(4) 

When compared to the best gradient deter- 
mined using M33 PNe (Eq. 3), we confirm that 
these gradients are consistent within the errors. 
We infer that the metalhcity gradient of M33 has 
not changed from the epoch of the formation of 
the PN progenitors to the present time. 

We can use for comparison the gradient of the 
best measured element, oxygen, or the weighed 
average of the slope of 0/H, Ne/H, and S/H gra- 
dient, and shown in Fig. [IH -0.033±0.003 dex 
kpc~^. 

In Table [8] we compare our metallicity gradient 
with the literature values derived using different 
stellar populations, a-element and iron-peak gra- 
dients are separated by a double horizontal line 
since their different origin (SNII vs. SNIa) sus- 
pects a different metallicity distribution through 
the disk and possibly a different gradient. In each 
group, a-elements or Fe/H, gradients are ordered 
by the increasing age of the corresponding stellar 
population. We considered only metallicity gradi- 
ents computed with a statically significant sample, 
thus, e.g., excluding the metallicity gradient with 
Cepheid stars (Beaulieu et al. 2006). 

It is remarkable that both the oxygen and neon 
metallicity gradients of M33 do not change within 
the errors even considering populations of different 
ages, though they are slightly different, the neon 
gradient having a slightly steeper slope. 

In the following we discuss which evolutionary 
scenario and chemical evolution models are ex- 
cluded or favored by our results. 

The first chemical evolution model attempting 
to the descri be the evolution of M33 was that by 
Diaz & Tosi (|l984 . Their model assumes a con- 
stant and uniform infall rate and fails to reproduce 
the metallicity gradient observed in M33. Other 
models that describe the origin and evolution of 
metallicity gradients are listed in Table 7, where 
we report the Reference (column 1), the model 
type (column 2), the resulting oxygen gradient 
(column 3), and its variation with time (column 
4), where a plus sign indicates a flattening of the 
metallicity gradient with time. 

The first three models listed in Table 7 are a 
family of multiphase models, differing mainly from 



the origin and rate of the infalling material. A 
multiphase model considers the galaxy divided in 
two or more zones, typically the halo and the disk, 
composed by baryonic material in different phases: 
diffuse gas, clouds, stars, and stellar remnants. 
These models evolve with time (see e.g. Ferrini et 
al. 1992, 1994 for a complete description), and are 
characterized by a flattening of the radial metal 
distribution with time, due to the inside-out for- 
mation of the disk. 

MoUa et al. (|l997l ) specifically analyzed the 



time-evolution of the metallicity gradient in sev- 
eral spiral galaxies, where the disk is formed 
by the primordial gas of the spheroidal proto- 
galaxy which collapses onto the galactic plane and 
forms out the disk. For M33 they obtained that 
A(0/H)/A(R)~ -0.21 dex kpc-\ flattening with 
time at a rate of -1-0.005 dex kpc~^ Gyr"-^. 

Molla & Diaz (jiooi) applied their results to 
M33 and found met allicit y gradients similar to 
those by Molla et al. ( 19971 ). Gradient evolution is 
not computed explicitly, but a clear flattening of 
the gradient is evident in their figures. 

M07b built the so-called "accretion" model re- 
producing the slow formation of the disk of M33 
from the inflow of intergalactic gas, which predicts 
A(0/H)/A(R)- -0.067 dex kpc-\ also flattening 
with time at a rate of +0.003 dex kpc~^ Gyr~^. 

Other types of models, such as those by Chi- 
appini et al. (1997, 2001), were not applied di- 
rectly to M33, but are notable since they predict 
steepening of the metallicity gradients with time. 
These models assume two main accretion episodes 
for the formation of the galaxy, the first form- 
ing the halo, the bulge, and the thick disk in a 
short timescale, and the second forming the thin 
disk, with a timescale that is an increasing func- 
tion of the galactocentric distance. Different as- 
sumptions on the star formation threshold and on 
the timescales of the infall episodes produce the 
gradients slopes and evolutions reported in Table 
[7] (see Chiappini et al. 2001 for a detailed de- 
scription). Note that the behavior of Model D is 
different from the other models by the same au- 
thors in that the assumed slower formation of the 
outer halo does not influence the disk formation, 
resulting in steeper gradients. 

If for example we assume that the age of all 
PNe is 5 Gyr, then according to the models by 
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Molla et al. we would expect a PN gradient which 
is -0.025 dex kpc~^ steeper than that of the H ii 
regions. Using the M07b models for the same PNe 
age, the PN gradient is -0.015 dex kpc~^ steeper 
than that of the H ii regions. On the other hand, 
with Models A, B, C by Chiappini et al. (2001) we 
would observe PN gradients that are flatter that 
those of H II regions by -h0.04, +0.02, +0.03 dex 
kpc~^ respectively. These variations could be even 
larger for older PN populations. In summary, only 
the mo dels b y M07b, and model B by Chiappini 
et al. ( 200lh . are marginally consistent with the 
observations. 

To match the small metallicity gradient vari- 
ation with time we observed with the models of 
M33 we could invoke a lower accretion rate from 
the ISM than that used in the M07b models. The 
resulting simulation should produce a flat metal- 
licity gradient, and a small slope variation in the 
last 8-10 Gyr, as observed. Detailed analysis of 
galactic chemical evolution models are out of the 
scope of this work. 



PNe in the disk of M33 are very similar to those 
of the LMC. The comparison between M33 PNe 
and H ii regions indicates a negligible global en- 
richment of the M33 disk from the epoch of the 
formation of the PN progenitors to the present 
time. 

The radial metallicity gradients of those ele- 
ments which are not modified during the lifetime 
of LIMS- oxygen, neon, sulphur- were derived. 
The best measured element, oxygen, has a slope 
-0.031 (± 0.013) dex kpc'^ in agreement within 
the errors with the same gradient derived from 
H II regions -0.032 (± 0.009) dex kpc'^. 

Since the metallicity gradients of PNe and H II 
regions are practically indiscernible and the mean 
abundances of the PN and H ii region populations 
are very close, we conclude that the chemical en- 
richment in M33 from the time of the formation of 
the PN progenitors to the present-time has been 
nearly negligible. 



8. Summary and conclusions 

We present spectroscopic observations of a large 
sample (102) of PNe in the spiral galaxy M33. The 
observations were secured with the spectrograph 
Hectospec on the 6.5 m telescope MMT. 

For 32 PNe the electron temperature was di- 
rectly computed. Empirical relationships among 
electron temperatures and bright emission lines as 
He II A4686, and [O ii] A3727 were used to derive 
electron temperatures of the remaining PNe. 

Abundance diagnostic diagrams were built to 
study the PN population. From the plot of N/0 
versus He/ H we found that the PNe that fulfill 
the Dopita (|l99ll ) definition of Type I PN are 19. 
Most PNe are non Type I, implying a population 
mainly composed by PNe from old progenitors, 
with <3M0 and ages >0.3 Gyr. 

A tight relationship between the 0/H and 
Ne/H abundances was found, excluding the mod- 
ification of both elements by PNe progenitors and 
ensuring oxygen to be a good tracer of the galaxy 
metallicity. 

The average chemical abundances of the PNe 
were compared to the PN population of the LMC, 
the SMC, and the Galaxy. Generally speaking, 
the elemental abundances of the a-elements of the 
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Table 7 

Chemical evolution models of M33 



Ref. 


Model 


A(0/H)/A(R) (present time) 


(A(0/H)/A(R))/A(t) 






(dex kpc"""") 


(dex kpc"^ Gyr"^) 


(1) 


(2) 


(3) 


(4) 


Molla et al. (1997) 


multiphase, iiifall from collapse of the halo 


-0.21 


-1-0.005 


MoUa & Diaz (2005) 


multiphase, infall from collapse of the halo 


-0.20 




M07b 


multiphase, accretion model 


-0.067 


-1-0.003 


Chiappini et al.(2001) 


model A 


-0.065 


-0.008 




model B 


-0.045 


-0.004 




model C 


-0.070 


-0.006 




model D 


-0.20 


-1-0.05 



Note. — (1) Reference; (2) main characteristics of the chemical evolution model: (3) slope of the present-time O/H radial gradient; 
(4) time variation of the O/H gradient. 



Table 8 

a sample of literature abundance gradients 



Element 


Gradient 


Population 


Age 


# 


Ref. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


O/H 


-0.05±0.025 


B giants 


few Myr 


30 


Urbaneja et al. (2005') 


O/H 


-0.032±0.009 


H II regions 


few Myr 


89 


RS08 + M07a 


O/H 


-0.039±0.033 


young PNc 


<0.3 Gyr 


19 


this paper 


O/H 


-0.031±0.013 


All PNe 


<10 Gyr 


91 


this paper 


O/H 


-0.030±0.013 


old PNc 


0.3-10 Gyr 


72 


this paper 


Nc/H 


-0.05±0.02 


H II regions 


few Myr 


30 


Willner & Nelson-Patel (2002) 


Nc/H 


-0.058±0.014 


H II regions 


few Myr 


25 


Rubin et al. (20081 


Nc/H 


-0.025±0.037 


young PNe 


<0.3 Gyr 


18 


this paper 


Nc/H 


-0.037±0.018 


All PNe 


<10 Gyr 


55 


this paper 


Nc/H 


-0.051±0.019 


old PNe 


0.3-10 Gyr 


37 


this paper 


Fe/H 


-0.06 


AGB 


~4-6 Gyr 




Cioni et al. (2008) 


Fe/H 


-0.04±0.02 


RGB 


>8 Gyr 




Kim et al. (2002) 


Fe/H 


-0.06±0.01 


RGB 


>8 Gyr 




Tiede et al. (2004) 



Note. — (1) Chemical element; (2) slope of the radial gradient with its error; (3) stellar population; 
(4) typical age of the stellar population; (5) size of the sample; (6) reference. 
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A. PN039: a planetary nebula with a [WC]-nucleus 



PN039 has a WR central star. Its spectrum shows prominent and wide features of He II at A 4686, 4541, 
5411 and of He I at A 4471, 5876, 6678, 7065. In addition other stellar Hues of [CHI] at A 4056 and 4649, 
and [CII] at 7107 and 7124 were detected. The FWHM of these lines is 15A, to be compared instrumental 
resolution of about 6 A. Thus, following the classification of Acker & Neiner ( 2003f) PN039 has a late- type 
[WC]-9 or [WCJ-10 nucleus due to the presence of very bright [CHI] and [CII] lines respect to the [CIV] lines. 

The ratio of the sulphur doublet indicates a low density limit <100 cm~'^, and [O ill], 4363/5007 ratio is 
completely anomalous for a normal PN, giving an electron temperature 30,000 K. The reddening of the 
nebula is null, within the errors, indicating a little presence of dust. 
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have improved the paper. We thank D. Fabricant for making Hectospec available to the community and the 
TAC for awarding us the observing time. We thank the Hectospec instrument team and MMT staff for their 
expert help in preparing and carrying out the Hectospec observing runs. We thank N. Caldwell, D. Ming 
and their team for the help during the data reduction. Thanks to Katia Cunha for enlightening scientific 
discussion, and Richard Shaw for his help with the nebular package. 
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